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1. Introduction 

High energy QCD has reached a mature stage of development, in its description of dilute-dense scattering. 
Deep inelastic scattering with nuclei provides a good example of this. The main physical phenomena 
that emerges for this type of scattering has been discussed, and the non-linear equations that govern such 
processes have been derived and discussed in detail [1-8]. On the other hand, the scattering of the dense 
system of partons with the dense system of partons has been actively studied [9-15], but with limited 
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success. This is in spite of the fact that this scattering is closely related to nucleus-nucleus scattering, 
which is the source of most of the experimental information, for the dense parton system. 

At the moment there exist two general approaches to high energy QCD: the BFKL Pomeron calculus 
[1,2,9,16,17], and the Colour Glass Condensate approach (CGC) [3,7], which lead to the same non-linear 
equations [5, 6] for dilute-dense scattering. The interrelation between these two approaches is not clear 
at the moment. However, the equations that describe nucleus- nucleus collisions have not been derived, in 
spite of considerable progress made in this direction [10,12,13], while in the BFKL Pomeron calculus, such 
equations have been proposed in Ref. [9]. These equations have been on the market for some time, but 
unfortunately, only three attempts to solve them are available in Refs. [18-20]. 

The main goal of this paper is to find the solution to these equations. In the next section we re-derive 
the equation of Ref. [9] in momentum representation, which turn out to be the most economical way of 
finding the solution. In section 3 we will find the semi-classical solutions to the equation, which describe 
dense-dense scattering outside of the saturation region. Section 4 is devoted to finding the solution inside 
of the saturation region. In the conclusion section we summarize our results. 

2. The BFKL Pomeron Calculus 

The goal of this section is to find the equation for nucleus-nucleus scattering in the momentum represen- 
tation based on the BFKL Pomeron calculus, based on the main idea of Ref. [9] that the equation for 
nucleus-nucleus collisions can be found from the equation of motion for Pomerons. 

2.1 Functional integral formulation of the BFKL Pomeron Calculus Sq 

The BFKL Pomeron calculus can be written through the functional integral [9] 

Z[$,$+] = j D^D^+e^ with S = Sq + Si + Se (2.1) 

where Sq describes free Pomerons, Si corresponds to their mutual interaction while Se relates to the 
interaction with the external sources (target and projectile). Here the free action is given by: 

Sq = j dY' j (fxi<fx2^^ {xi,X2,Y')VlVl[-^ + H]^{xi,X2,Y') (2.2) 

Define the following Fourier transform, 

$t(xi,x2,y') = $+(xi2,6,y') = xl2 j Skie-'^^-''^^^^ {ki,b,Y') (2.3) 
$(xi,X2,y') = $(a;i2,6,y') = XI2 J d\2e"''''''^ {k2,b,Y') (2.4) 

Consider the free action as the sum of two terms: 
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•5*0 — 'S'o + ^0 (2-5) 

where: 

So' = J dY' J d^xi(fx2^Hxi,X2,Y')vM-^^{xi,X2,Y') (2.6) 



This can be re- written as: 



S^ = 4 1 dY' I d%d'xu I d^fcid^A^ae^-^i^c&t {ki,b,Y ') xj.vM (^?2 e"''^"^^^^ (^2, 6, 1^') ) (2-7) 
where in the last step the replacement d'^xi <f'X2 = A<f'hd'^xi2, where, 



= ^ Xi2 = Xi-X2 (2i 

where b is the the impact parameter. In two-dimensional polar coordinates: 



2 



Let us introduce the next change of variable I = In fc^ , such that the two-dimensional Laplacian V^^ 
with respect to /c2, simplifies to the following differential operator: 



Vfe = 4e-''=.|^ + e-''=.^ (2.10) 
According to the definition of Eq. (|2.10| ), the next term icf2V^^ V^^ (^12 e^*'^^'^^^) can be recast as, 

2 ^2 ^2 ( 2 -ik2-xi2\ — ^2 ^4 -ik2-xi2 

^12 \^a;i2 e j — Xi2^ j^^w ^^^e 

— --r^ ( i.4 p-«A:2-a;i2 

— ^12Vfc2 \^'''2 6 

— Vfc2 P2 Vfc2 e 



9Z / 86 \dP de 



16 I 1^ + 1) +7^7:1(^ + ^1 (2.11 
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where V^^ is the two-dimensional Laplacian derivative with respect to ^2- Inserting Eq. ( ^.11| ) back 



into Eq. (2/7), leads to the expression: 



5o = 4 j dY'j (fbcfxuj (fki(fk2 (2.12) 
X e*'=---ci>t (fci, 6, y ') { V^, (ki Vl e-^^--- ) } {k^, b, Y') 



32 / dV / d%d^xi2 / d^kidOdle^ (2.13) 



where in the last step, the new variable / was introduced. We use the convention that derivatives only 
act on terms inside of the curly brackets {}. Hence in Eq. ( 2.13D derivatives only act on e~**'2-2:i2_ Assuming 



that the ^-dependent part of the integrand is a purely periodic function of 6, then the derivatives can be 
re-ordered through integration by parts, to yield: 

Sq = ?,2 j dY' j d%d^xi2j d^kidOdle^ e'^^-''^^ e-'^^-""^^ (2.14) 

x*.(^.M^')(ff|.i)%J^)(^.|)^*(^..nl 



32 



j dY'J dHd^xuJ d^kidOdle^ e'^'-'"'^ e-*fc2-xi2 ^2.15) 

.,.(.,,v'){(|+i)'^|,.(.,..,n} 

where in the last step it was assumed that <I> {k2,b,Y') is purely a function of = and not a function 
of the angular coordinate 9, such that ^-derivatives vanish. Now returning to the integration variables 
J dOdl = 2 J" d^k2 where I = In/c^, and integrating over xi2 leads to the delta function (27r)^ 5^ {ki — k2). 
The delta function is absorbed by the integral over A;2-space resulting in the expression: 

5o = 64(2^)2|dy'y"d2^y"d2^i<I>t(A;i,6,y')|(| + l)' (2-16) 

where I = In/c^ This part of the action gives the following contribution to the equation of motion, 
which stems from the condition 5Sq/ 5^'^ {k,b,Y) = 0: 
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6S^/S^Hk,b,Y) = QA{27rf(^^^+l^ ^±^{k,b,Y') (2.17) 

Recall that in the above calculation, for the sake of simphcity we considered 5*0 = Sq+Sq, and we calcu- 
lated the Fourier transform of the Sq part. The full variation 6So/6^^ {k, b,Y) = 6 {Sq + Sq) /6^^ {k, b, Y) 
is given by the expression: 



SSo/d'^HkAY) = 64(2.)'(| + l)' ^ (^-^)^(^'&'^') (2-18) 

The general properties of the Hamiltonian H have been discussed in Refs. [9, 17] in coordinate repre- 
sentation, which read as follows in momentum space representation: 

Hf{k,Y) = ^ J d^lKikJ) {fik,l,Y) - fik,Y)} (2.19) 

K{k,l) = ^^-^ (2.20) 

^ ' ' ll(k-l)l ^ ' 

2.2 Interaction term of the action Sj and the Field equation 

The mutual interaction of the Pomerons is described by Sj. The interaction which is related to the external 
sources (target and projectile) are not consider in this paper. In this approach Si is given by: 



Sj ='-^fdY' [ ^ {{L,,<^ {x„x„Y')) c|>t (,„,3,y-) ^,,,,,,Y') (2.21) 



27ra| r ^y, f (Pxid?X2d?xz 

Nc J J ^12-^23-^31 

+ (Li2$^ (xi,X2,y')) ^ (^2,X3,y') $ (X3,xi,y')} 

= '-^ [dY'f ''yi?lf''' {{Li2^ (xi,X2,y')) (-2,X3,r') (X3,x„y') (2.22) 

JVc J J •''12-^23'''31 

+ (Li2$^ {xi,X2,Y')^ ^ {X2,X3,Y') $ (X3,xi,y')} 

where in the last step the replacement d'^xid'^X2 = 4d^6(i^xi2, (see Eq. ( p.8D ), and where the following 
differential operator was introduced: 

L,2 = xt2'^lVl (2.23) 
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Now the above defined Fourier transform of Eqs. 
its assumed that b » xi2, where b = {xi + X2) /2, i.e. 
of the dipole. Inserting the Fourier transforms of Eqs. 



(U) and (U) are inserted into Eq. (|]22|), where 
the impact parameter is much bigger than the size 
(|]|) and into Eq. (|^) gives, 



2Trai 



Si = —-^ / dY' / 'i(fb(fxi2(fx3 / (fki(fk2(fk 



exp {ik2 ■ X23 + iks ■ X31) 



(2.24) 



X 



12 



+ exp {-ik2 ■ X23 - iks ■ X31) 



'L12 (x?2 e^'^i-^i^y 



^12 



$ {ki,b,Y') $t (^2, 6, y ') $t (^3^ i,^Y') 



cDt (A:i,6,y')cI>(A:2,6,y')^(A:3,6,y') 



The integration over X3 leads to the Dirac delta function (27r)^ (5^ {k2 — k^), where 5^ {k2 — k^) labels 
the delta function in two dimensions. The delta function is absorbed by the integration over k^ which leads 
to: 



Sj = ^^ (27r)2 / dV I 4dHd^xi2 I Skid^k2 
I exp {-ik2 ■ X12) 

+ exp {ik2 ■ X12) 



(2.25) 



$(A;i,6,y')cl>t {k2AY')^^ {k2,b,Y') 

•^12 / 

^^^^^^^^^ '^HhAY')^{k2,b,Y')^ [k2 , 6, y ) } 



According to the definition of Eq. ( |2.23| ), the term L12 {x\2 exp {—iki ■ X12)) / x\2 can be recast as, 



Ll2{ 



„2 ifci-xi2j 



Xi2^ 



^12 



2 ( 1,4 V72 ^-jA:i-xi2 



2 V72 V72 / 2 -ifci-xi2 
■^12 ^a;2 \ -^12 ^ 



V^, (A;? Vi^ e 



9/2 + 06' 



-ifcl-xi2 



(2.26) 



(2.27) 



where I = InA;^ and where the last expression was discussed and introduced in the first section (see 
Eq. ( ^.llD ). Inserting this result into the interaction action of Eq. ( 2.25 ), after some algebra one arrives at 
the expression: 
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Si = 16 (2^)2 / dY' I A(fb(fxi2 I (fki(fk2 (2.28) 



exp (-i (A;i + A:2) • X12) (^ + (^i, 6, F') ^^^2, 6, 1^') ^U^2, ^ i^') 

+ exp (i (fci + fcs) • X12) + 1) {kiAY')^{k2AY')^{k2AY')^ 

where it was assumed that the functions $ (fci, b,Y') and (/ci, 6, y' ) are purely functions of the radial 
k\ coordinate and do not depend on the angular coordinate 9, such that derivatives that appeared in 



Eq. (|23^ ) have been dropped. The integration over X12 leads to the Dirac delta function (27r) 5"^ {ki + h 



2) 



The delta function is absorbed by the integration over k2 which leads to: 



Sj = 16 (27r)^ I dV I Acfb I d^ki (2.29) 



l^<^^-k„b,Y')^H-k„b,Y') (l + l) ^^{k^AY') 

<^ {-k,AY') $ {-k,AY') (I + 1)' ^2^^ [kiAY')^ 



where I = Infcf. Integrating the last term in Eq. ( 2.29| ) by parts, and taking into account that 
d^ki = ^dOdle^ where 9 is the azimuthal angle, Eq. ( 2.29| ) simplifies to the following expression: 



Sj = 16 (27r)4 I dV I Ad^b I d^ki (2.30) 



I ci>t (-k„b,Y') ci>t {-k„b,Y') (I + 1)' ih,b,Y') 
The last equation Eq. ( p. 301) allows us to study the classical equation of the action, which can be 

ss §s 

obtained from the conditions [9,201: - — -; — ; = and ^ , — ; = 

6^{k,b,Y) 6^T{k,b,Y) 

2.3 Classical equations 

The functional derivative of the action for the effective pomeron field theory, with respect to {k,b,Y) 
can be found from the results of Eqs. ( p. 18 ) and ( 2.3C| ), which give: 
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where / = In A;^. The equation for nucleus- nucleus scattering can be derived from the following equation 
of motion [9]: 



5{So + Si) /5^^k,b,Y) = and 5 {Sq + Si) {k,b,Y) = (2.32) 
Now the following approach is used to average these equations: 

f D<^>D^Wik,Y;b) e^i'^''^') , , 

0{k,Y;b)) = ^ ^ ' ' \, (2.33 

Introducing the following new functions: 



N{k,Y;b) = 2TT^as(^^{k,Y;b)j N^k,Y;b) = 2TT^as(^<^Hk,Y;b)J (2.34) 

then the equation of motion of Eq. ( 2.31| ) reduces to: 



+ asl2NU-k,.b.Y') + l)' ^ N(k.b.Y') + (^l + l)' ^ N^-k.b,Y ') \ 

Using the following property of the BFKL Pomeron, namely, 

we can derive the second equation of motion, by taking the functional derivative with respect to 
This generates the following equation of motion: 



dl ^ dP \ dY 



as^N i-h,h,Y') 



d_ 

di 



NHki,b,Y) (2.37) 



In Eq. (p.35|) and Eq. (2.37) the following identities were implemented: 



2 



^'ik,Y;b)) = ^{^'^ik,Y;b))) (2.38) 
$ {k, Y; b) (A;, Y;b)) = {k, Y'h)) y. {k, Y; b) 



For a discussion on why these equations are correct, within an accuracy of about in the case of 

nucleus-nucleus scattering, see Refs. [9,20]. 

3. Semiclassical solution 

3.1 The system of equations: general approach 



In this section we find the semiclassical solutions to Eq. ( 2.35| ) and Eq. ( ^.37 ). In the semi-classical 



approximation, we are searching for solutions of the following form [21]: 

iV(/ = ln(A;2),6,y') = exp (S (/, 6, F')) where S = ujY' - {I - 7) / + /3(6) 

iVt(Z = ln(fc2),6,y') = exp(5^ (Z,6,y')) where = J Y' - {I - + /3\b) (3.1) 

where uj{l, b, Y') and ^{1, b, Y') are smooth functions of Y' and I, and the following conditions are assumed: 

Mipri « Mhhn « .(,,6,r)(i-7ai.,i")); (3.2) 

'-^^ « (l-7«,«"))^ « <^iI,i.Y'){l-,iLi,Y')): (3.3) 

with analogous conditions for the functions io^ {l,b,Y') and 'y^ {l,b,Y'). Assuming that for Eq. (p.l|), the 
method of characteristics can be applied (see, for example, Ref. [22]) to solve the non-linear equation. 
Notice that 

He^ = a5y"i^(fc,A;')e^(W) (34) 
where xil) = 2'0(1) — (7) — ^"(1 — t) (3-5) 
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where ip{z) = dliiT{z)/dz is the Di-Gamma function. In the semiclassical approach (for 7 a smooth 
function), inserting the definition of Eq. ( |3.lD into the equations of motion of Eq. ( p.35 ) and Eq. ( p. 37 ), 
and using the conditions of Eqs.(^) - (|3.3|), leads to the following formulae in the semi-classical approach: 



w - asx (7) + ^ [l^j e'^^ + ase^ =0 

-uj^ - asx[l^^ + i^il) Ois e^'' + ase^ = (3.6) 
where the following functions were introduced with definitions: 

""^^^ = 4(27-1)2 ^ = S-luKij) 5t = St-lnK(7t) (3.7) 

For the equation in the form 

F(y',/,S,7,w) = (3.8) 

where S is given by Eq. (|3.1|) , we can introduce the set of characteristic lines on which l{t),Y'{t), S{t), oj{t), 
and 7(t) are functions of the variable t (which we call artificial time), that satisfy the following equations: 

^'•^ d-t= =-"^^ 
dY' 

(3-) -77 = 7^7 + ^Fu. = as (1 - 7) + w 

(xt 07 

(4.) ^ = -(F^ + 7i^s) = a5^(7+) (l-7^)e^' + a5(l-7)e^ 

(5.) ^ = -[Fy + u:Fs) = -asw^K (7+) e^^ - a^we^ (3.9) 

where F; = ^ — '^'^ — ^, and the terms in the first of Eqns. (^) that depend on S*^, (^"^ 7^, are treated 

as if they depend explicitly on Y' and /. The same five equations we can write for the second of Eqns. 
which have the following form: 





lit 


~- = d7t 


(2.) 


dt 


= ^..t = -1, 


(3.) 


dS^ 
dt 


= 7^^^^ + ^^Pu.^ = 


(4.) 


d-f^ 
dt 




(5.) 


duj^ 
dt 


= -{Fy + UJ^ Fs) 



as (l 



■asu;K{j)e^ -asoo^e^^ (3.10) 
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3.2 The system of equations: linear approximation 



We start with the solution of Eq. ( ^.351 ) and Eq. ( 2.37| ) in the kinematic region where the non-linear 



corrections are small. In this region, we can reduce this system of equations to the solution of the BFKL 
equation, both for b, Y') and A^^ (/, b, Y'). Actually in order to find the solution, it isn't necessary to 
solve the equation, since the t-channel unitarity constraints for the BFKL Pomeron can be used instead 
(see Refs. [1,16,23]). This is given by: 

NBFKL{L,b,Y) = J (fb' J dlNlpj^^(^L-l,b-b',Y -Y'^ NBFKL{l,b',Y') (3.11) 
Using the explicit form for the BFKL solution [16, 17], then Eq. ( 3.1l| ) reduces to the following expres- 



sion: 



NBFKL{L,b,Y) = y ^nBFKL(7,ft) (3.12) 

where L = ln(A;p/A;|) and / = ln(/c^/fe|). kp and kt are the momenta of the dipoles in the projectile 
and the target, respectively. The function hbfkl (7, ^) is determined by the initial condition at y = 0. 

The integration over I leads to the delta function (27r) 5 (7 — 7')^ ^ind hence the 7' integral is imme- 
diately solvable by setting 7 = 7' everywhere in the integrand. Integrating over the impact parameter b' 
yields: 

j d%' n^BFKLil^ b - b') n^^^^(7, b') = riBFKdl. h) (3.13) 
In light of this Eq. ( ^.12 ) simplifies to 



NBFKL{L,bX) = I e^^^^'^^''-^'-'^^'^nBFKLh,b) (3.14) 



Incidentally Eq. ( ^.14 ) satisfies Eq. ( 3.11J ). Therefore, we have reproduced Eq. ( |3.11| ), where and 



are given explicitly by: 

NU,{L-l,b-b',Y-Y') = I ^^ni,^,h,b-b')e-sxMiy-y')-i^-.)iL-i) (3.15) 

NBFKLil,b',Y') = [pl,nBFKLil',b')e^^^^^'^^' + ^'-^'^' (3.16) 
J 2tti 

where the functions n^BFKii'y-' b~ b') and ubfkl {i\ b') can be found from the initial conditions for N^bfkl 
atY' = Y and Nbfkl at Y' = 0. 
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N(L,Y) 




Figure 1: The BFKL Pomeron: i-channel uni 
tarity. 



The integral expressions of Eqs. ( 3.15 ) and ( ^.16 ) can 
N\ L-I,Y- Y') be used to derive the exphcit expressions for 

^BFKL (^L-l,b-b',Y - and Nbfkl {I, b', Y') in the 
semi-classical approach, by integrating using the method of 
steepest decent. This is assuming that n^BFRL^'y^ h — b') and 
nsFKL {'y' ,b') are continuous, and differentiable functions 
of 7 and 7'. 



N(I,Y') 



The saddle point equations, for the functions in the exponent in Eqs. ( ^.15 ) and ( 3.16 ) are: 



^^dxpspl^y_y.^^^_l ^ oanda^ 



d-fsp " " - o ^^^^ 

The integrals of Eqns. (3.15) and ( 3.16| ) can be solved using the method of steepest descents to yield: 



(3.17) 



Nlpj,^{L-lb-b\Y -Y' 



(3.18) 



2-K 



asx"{lSp){Y-Y') ^^^^^ 
= nBFKL{lSP,b-b')e^^ 



{7sp,b-b') exp (asxhsp) {Y - Y') - (L - /) (1 - jsp)) 



Nbfkl [hb-b'X 



(3.19) 



27r 



^sx"{l'sp)Y 
= flBFKL [I'sp^b') 



- nBFKdl'sp, b - b') exp [asx [I'sp) Y'p-l{l- 7^p)) 



S = as{x{l') - i^-ysp)'-^)y'-^ = ^s{x{lSP)-{l-lSP)^^^) {Y-Y') 

where all slowly changing terms, have been absorbed by the functions and h. Eq. ( [3.9| ) together with 
Eq. ( |2.35| ) have the following form in the linear approximation: 

(l-)^ = -a5-^; {2.)uj = asxiin (^"^ ^ = ~ ^) + ^! (4-) 3^ = 0; (3.20) 



dY' 



It is easy to see that Eq. ( |3.20| ) leads to the same S as Eq. ( p. 18 ) and Eq. ( p. 19 ). One can see that for 
7 = 7cr for which 



X{lcr) + (1 - Icr) = 



(3.21) 
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L-aY 



Figure 2: The trajectories of the Imear equations for and N: L = Inyk^/kfj ,1 = hi(fc^/fc^) and a = 

OisX ilcr) her- Red lines denote the critical trajectories. Only the trajectories to the left (for TV^ ) and to the right 
(for TV) are shown. 

Then 5 = 0. The equation for this line takes the form: 

I = as^^^Y' (3.22) 

1 - 7cr 

Repeating ah the steps of the calculations above for S"^ , we find that 5^ = on the line 

L-l = as^P^ {Y-Y') (3.23) 



using the same jcr from Eq. ( 3.21 ). The general pattern of trajectories for the linear equation is shown in 
Fig. |2[ It should be stressed that for Eq. (|3.9D and Eq. (^.lOl) , we have different trajectories but they are 
parallel and shifted by Al = as (x (7cr) /jcr) Y. 

We need to know the initial conditions for S{S^) and 7 (7''^) at y' = {Y' = Y). We choose the 
McLerran-Venugopalan formula [3] which is written for the dipole-target amplitude and is given in terms 
of N as: 

N{r,b,Y' = 0) = 1 - e^p{-r^kf{b)/4) (3.24) 

where the initial characteristic momentum kf{b) oc Ta (6), and 7a (6) is used to denote the number of 
nucleons inside the nucleus with fixed impact parameter b. In momentum representation N is equal to (see 
Fig. §: 

No{k,b,Y' = 0) = y" rdr Jo (A;r)iV(r, 6, y = 0) /r2 = lr(0,r=^) or 5 = In (Iro(r)) 



d\n(N{k,b,Y' = 0) 

Fq (t) that appears in Eq. ( |3.25| ) is the Euler incomplete gamma function ( see formulae 8.50 in 
Ref. [24]). 
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Caution should be taken here, since we cannot trust the McLerran-Venugopalan formula at small 
dipole sizes. Indeed, we know that in the limit of perturbative QCD, 70 is 7 — ?> contradicts this formula. 



The problem is that the simplified version of Eq. (3.24) does not reproduce the perturbative QCD limit 

^ at r ^ 0, while the correct behaviour should be 



at r 



oc r 



^ 0. The relation in Eq. leads to N 

N (X r"^ Inr^. Nevertheless we will use Eq. (3.24), because our main interest lyes in the region in the vicinity 



N(T), 7(T) 




of the saturation scale, where Eq. (3.24) reproduces the amplitude for N quite well. 

For N'^ the initial conditions look the same, but = 
P /kj. Inotherwovds, N'f {k,b,Y' = Y) = e^o = iro(rt). 
One can see that using Eq. (^) (4) and Eq. ( |3.10| )(4), we 
obtain that 5t(Z,y') = S {L,Y) - S{l,Y') + Sl It 
means that uj^ = —oj and 1 — 7"'" = — (1 — 7). It is easy to 
see that the system of equations in Eqns. ( ^.lOD degenerate 



to the system of equations in (|3.9| ). 



3.3 The final system of equations 



Figure 3: Initial conditions: 7Vo(r) (solid line) Therefore, instead of Eq. (pj) and Eq. (H) , we can solve 
and 7o (dotted line). 1 p n • 

the followmg system of equations. 



(3.) 

(4.) 
(5.) 



^ = «5 (1 



dl 
It 
dS 

lit 
doj 

Ht 



as ■ 



dj 



(2.) 



dY' 

IT 



1, 



-asK{^) (l-7)e^(^'^)-^(''^') + 5o + as{l- 7)e^(''^'), 



(3.26) 



For the sake of simplicity, the label has been omitted everywhere in Eq. ( ^.26 ). Eq. ( |3.260 can be 
re-written in a different form, namely 



dl 

dY' 



as 



(2.) 



dS 



as (1 



d-f ' ' dY' 

^SiL,Y) - S(l,Y')+So + eS(l,Y) 



-,)^^+.. (4,) p. 



1-7 



(3.27) 



Eq. ( |3.27| )-4 has the solution w 
conditions. 



dx{l) /dl + a;/(l - 7) 

as Const{\ — 7) where Const can be determined from the initial 

x(7o)-«(7o) e^^^'^^ - e^o 



Const 



1-70 



(3.28) 



Notice that the value of this constant depends on S{L,Y), but for large and negative S{L,Y) which 
we are dealing with, for 70 < 7cr, we can safely neglect this dependence. Introducing S = S {l,Y') — 
{S (L, Y) + 5*0) / 2 we can rewrite Eq. ( \i.27 )-5 in the form 
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Fig. |-a Fig. |-b Fig. ||-c 

Figure 4: Semiclassical solution: trajectories I (Y) (Fig. ^a), 7 versus asY (Fig. ) and S versus 0:51" (Fig. ^c). 
The solution with 7 close to jcr is shown in red by dotted line. asY is chosen to be equal 6. 



dS 



,iS(L,Y)+So)/2 



-k(7) 



■ S(l,Y') 



+ e 



dxil) /d'y + Const 



(3.29) 




We follow the following strategy for solving this equation. 
First we assume that we are looking for the solution in the interval 
{So, Smax}, where we have selected Smax- Using the initial con- 
ditions, we can calculate S {l,Y = 0) through and S^axj 
specify the coefficient in front of Eq. ( |3.29 ). Solving Eq. ( ^.27 )-2 
we will find Y' as a function of S. In particular Y' (Smax) 7^ Y, 
however trying several times we find a value for Smax that will 
give Y' (Smax)- Fig. ^ shows the solutions for four values of 70 
that are smaller than 'jcr- One can see that the trajectories and 
from Eq. (|3.3l|) (thick line) and the values of S and 7 are close to the solution of the linear equa- 
0.37 (solution of Eq. ( ^.21 ), thin tion. Indeed, dS/dY' < for 7 < 7^ and therefore, the value of 

S decreases due to evolution. Starting with 5o < one can see 
that the contribution of the non-linear terms becomes less and 
less at higher values of Y' than at Y' = 0. Therefore, only at small values of Y' we can see a deviation 
from the solution to the linear equation, as shown in Fig. ^ (see Fig. ^b for example). Eq. ( t3.27D -4 has 
the form 



Figure 5: Dependence of S {Y' = 0) = 5*0 

on 7, 

7cr = 

line). 



d'y 
dY' 



as (1 - 7) e 



{S{L,Y)+So)/2 



S{l,Y') 



+ e 



S{l,Y- 



(3.30) 



and due to the smallness of the factor exp {{S {L, Y) + ^o) /2), then d'y/dY' turns out to be small at large 
Y', leading to a constant 7, as shown in Fig. ^b. The value of 7cr is slightly different from the one from 
Eq. ( 3.2l| ), due to a contribution from the non-linear terms to Eq. ( |3.27 )-3, and it depends on the initial 
condition for S (Y' = 0) = Sq. Indeed, for the trajectory on which S {l,Y') = Sq is constant we have the 
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following equation: 



(l-7cr)^^7^ + X(7cr) " {^{icr) + 1) 6^° = (3.31) 
«7cr 



The dependence of Sq on from Eq. ( |3.3ip is shown in Fig. ^. One can see that for 5o < [Nq < 1), 
the value of •jcr is close to the solution to Eq. ( ^.21 ). From Fig. ^ one can see that we start the evolution 



from the value of 7 that is close to 7cr, 7 steeply increases to 7 = 7cr and freezes at this value leading to 
constant S almost in the entire kinematic region of asY'. For 7 > ^cr, the solution will lead to S {l,Y') 
that increases with Y', and we need to search for a different method of finding the solution, other than the 
semi-classical approach. 



4. Solution inside the saturation domain 



4.1 General consideration 



As discussed above, the semi-classical approach cannot be used inside of the saturation region. It should 
be mentioned that at large /, the amplitude behaves as ^l, which is certainly not the function for which 
we can use the semi-classical approach. However, it has been noticed in Ref. [21] that introducing new 
functions : (p (^'1 ^'1 b) and 0^ (/', Y', b) instead of N and A^^, defined as 



N 



{l,Y',b) =\ldl'{l- e-*('''^'''')), N^l,Y',b) = 



-</.t(r,y',6) 



(4.1) 



Then we can indeed use the semi-classical approach for these functions. The first observation is that, 



from the property of Eq. (3.4) and the definition of Eq. (3.16): 



^ dN{l^Y',b) ^ l^^^^^^^)^^^^^_^^^Y',b)e^p{asx{l)Y'-il-^)l) 

= j ^.[x (7) - (7 - 1) n (7, Y', b) exp {asx (7) Y' - {I - 7) /) - asN (/, Y' , b) 

f d\ dN(l,Y',b) ^ . , 

" ai) m - 

The second observation is related to the function C (7) , namely that its expansion with respect to 
(1 — 7) starts with (1 — 7)^ as the first non-zero term. That is: 

.f_d\ _ d^jz) , _ 1 dV(z) , _ 

^[ dl ~ dz^ ^'='dP 12 dz^ '"='9/4 ^^-^^ 
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where il^{z) = dlnT{z)/dz is the Euler ^/^-function (see Ref. [24]). Following the definition of Eq. ( [4.1D and 
assuming that cf) and (j)^ are smooth functions, we can replace * 

n-l 



2 \dl 
1 f d 



^YNUl,Y',b) 
dl ^ ' ' ^ 2\dl 



1 

2 V 5/ 



(4.3) 



71-1 



1 - e 



1 



n— 1 



5Z / 



(4.4) 



Inserting Eqs. ([4. 2)) , (|4.3|) and ( [4.4]) and into Eq. ( ^.35 ) and Eq. ( 2.37 ), and introducing the notation 
d(j)/dY' = uj^ and d(j^ /dY' = uj\ we obtain: 



= {io^+asC (-7^) ) (1 - 7<^)2 7<^e-^(^'^'^' ) + 2 (| + l) | ( e-^^''"'^' ) (/, 6, ) ) (4.5) 
= + «5/:(-7i)) (l-7i)%ie-'^(''^'"') + 2a.(| + l)' |{e-*^('A.' ) iVt (/, 6, y' ) } (4.6) 



-9it(i,;,,y') 



One can see that 



9/ _ 



l,b,Y' ) 



iV 



{l,bX)} = e-*{i (1- 67,^ + 77|) e-^ - i(l -47,^ + 37104.7) 



7,^(1 -7<A)'^(^^i"')} 



Dividing both sides of Eq. (4^) by age and introducing the new variable uj = oj/as, which 
corresponds to the change Y' to Y' = asY' we obtain: 



-O^ - C (-7,^) + 2 TVt (z, 5, y') _ V / _ ^ + 2iV (/, b,Y') =0 (4.5 

A similar equation can be written for cf)'^ . However, we assume that cf)^ (/, b,Y' ) = cf) (L — l,b,Y — Y' ) 
based on our experience with the semi-classical solution. Looking for the solution that has a geometric 
scaling behaviour [25], we expect that 4>{l,b,Y') is a function of one variable: 



In (Q2 (y',5) /k^ 



1 - 7cr 



(4.9) 



*For d4>/dl {d(j)'' /dl) we use the notations 7,/, and 7^. We hope that it wiU not lead to any confusions with the notation, 
due to the similarity with 7 and 7^, that was heavily used in the previous section. 
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4.2 Asymptotic solution 



Plugging in cj^ (l, b,Y' ) = (f) {L — l,b,Y — Y' ) and the geometric scaling behavior of Eq. ( [4. 8] ) that we 
anticipate, and assuming that (j) ^ 1 we obtain: 



G(7) = ZY 
where G (7) 



(4.10) 



Xilcr) ~ ^ ^/-^^ 1+47 + 37' 

7 + £ (7) + 



ZY 



1 - 7cr 

In (Q2 (y)/fe 



(1 + 7)27 



Xilc 



Y 



L 



It should be noticed that 7 is defined as 7 = 
see that we have three types of solutions to Eq. ( |4.10D : 



1 - 7cr 

z. The function G (7,^) is shown in Fig. 



One can 
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Figure 6: Function G (7) (solid line) versus 7. The red line shows the r.h.s. of Eq. ( [1.10| ) at zy 
shows G{'y) as it is written in Eq. ( 4.10 ), while 6(7) with £ = is plotted in Fig. ||-b. 



20. Fig. B-a 



1. At large 7 then G (7) — > xilci") / — lcr)l and we have the solution: 7 = (1 — ^cr) /xilcr) zy 
which translates into (j) = ((1 — Jcr)/x ilcr)) zy z; 

2. We have solutions at 7 — >■ n where n = 1,2,... which lead to = nz. We only need to take 
into account n = 1, since other values of n give smaller contributions at large z. In vicinity 7 — > 1 
G (7) = -1/(1 - 7) and the solution to Eq. { \1.10{ ) gives </> = (1 + l/zy) z; 

3. Solutions where 7 — )• —n we do not consider, since they lead to decreasing cp at large Y. 

It is interesting to notice, that at in the case where we restrict ourselves to the leading twist contribu- 
tions to the BFKL kernel [26], only the first solution survives (see Fig. ^b). 

Using Eq. ( [4.lD we can obtain the asymptotic solution for the scattering amplitude, namely 



N{zy) 
N{zy) 



1 ( 

2 io '^^^ ~ ~ 7cr)/X (7cr)) ZY z) 



ZY 



dz 1 



exp I 



i — leading twist; 
general asymptotic behaviour; 



(4.11) 
(4.12) 
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It should be noticed that the expression in parentheses acxtually gives the scattering amp-htude 
in the coordinate representation. 



4.3 General equations 

In general, Eq. (4.8) can be reduced to the differential equation by taking derivatives with respect to z on 



both sides of the equation. The resulting expression is: 

^5 = e-*W - e-<^(--), ^=;^, (4.13) 
dz dz 

This equation belongs to the class of delay differential equations, and the solution to this equation we 
hope to publish in an upcoming paper, since it is a separate and rather difficult problem beyond the scope 
of this paper (see for example Ref. [27]). 



5. Conclusions 



In conclusion, we summarize our results as follows. First we re-wrote the action of the BFKL Pomeron 
calculus, and we derived the equations in momentum representation. It turns out that the equations that 
we obtain have a simpler form, than in coordinate representation in the format that they were originally 
derived in Ref. [9]. Second, we found the semi-classical solution to these equations, to the right of the critical 
line . In the saturation domain, we reduced these equations to the class of delay differential equations, and 
we found their asymptotic solution. This solution shows, that the nucleus-nucleus amplitude in coordinate 
representation approaches unity as 1 — AA^, where AA^ oc exp (— zy), where zy = In (Q^ (Y) rf) , where 
is the size of the colourless dipole in the nucleus. If we take into account only the leading twist part of 
the BFKL kernel, the behaviour of AA^ is similar to the behaviour of the amplitude of the dilute-dense 
parton system interaction, given in Ref. [26]. 

We hope that this paper will show, that the problem of the nucleus-nucleus interaction in the framework 
of the BFKL Pomeron calculus, can be solved. We hope that this development will motivate further efforts 
towards understanding the dense-dense scattering system. 
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